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^D ' Abstract. We discuss calculation of nuclear corrections to the structure functions for the deep- 

Cn , inelastic scattering of muon and (anti)neutrino. Our approach includes a QCD description of the 

^ ■ nucleon structure functions as well as the treatment of Fermi motion and nuclear binding, off-shell 

Q I correction to bound nucleon structure functions, nuclear pion excess and nuclear shadowing. We 

^^ . emphasize the dependence of nuclear effects on the type and C-parity of (anti)neutrino structure 

functions. We also examine the interplay between different nuclear effects in the Adler and the 

Gross-Llewellyn-Smith sum rules for nuclei. 
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^ ■ Introduction 

^, 

The study of nuclear effects in deep-inelastic structure functions (SF) has by now 
a long and rich history. The measurement of the ratio M2{A/B) = F2/F2 of SF of 
two nuclei (usually a complex nucleus to deuterium) in DIS experiments with charged 
i7^ I leptons (CL) show that significant nuclear effects are present even in the scaling regime 

GN ■ ("EMC effect") [1]. Many theoretical and experimental studies of nuclear effects on 

structure functions followed, and an extensive phenomenology has been developed, even 
though many questions about the origin of the effect still remain. Most of the direct 
[^ I experimental information on nuclear effects in the structure functions come from the CL 

O ■ scattering data. Future measurements with neutrino (v) and antineutrino (v) beams are 

particularly interesting and would provide complementary information to CL data [2, 3]. 
The quantitative understanding of nuclear effects in neutrino scattering is of primary 
importance for a correct interpretation of experimental results and for the evaluation 
of the corresponding uncertainties. An accurate account of nuclear effects is impor- 
tant in the determination of the parton distributions and the high-twist terms [4, 5], the 
electroweak parameters in neutrino scattering experiments [6], and also for the under- 
standing of neutrino masses and mixing. The next generation experiments would im- 
ply precision measurements to disentangle small oscillation signals from neutrino and 
antineutrino interactions on nuclei. This would require to improve our knowledge of 
v(v)-nucleus cross sections in order to reduce systematic uncertainties. 

In this communication we summarize the results of recent studies of nuclear DIS with 
charged leptons [7] and neutrino [8], which were aimed to develop a consistent and 
quantitative model of nuclear DIS applicable in the analyses of existing data and in the 
interpretation of future experiments. 



Nuclear structure functions in charged-lepton DIS 

The approach of Ref.[7] involves the treatment of a number of mechanisms in the 
regions of large and small values of the Bjorken variable x. The physics scale of the 
separation between "small" and "large" x regions comes from the comparison of a 
characteristic longitudinal DIS correlation length \/{Mx) in the laboratory system (see, 
e.g., [9, 10]) with an average distance between bound nucleons. The DIS time is small at 
large x {x> 0.\), that justifies the approximation of nuclear DIS in terms of incoherent 
scattering off bound protons and neutrons (the impulse approximation, or lA). Important 
nuclear corrections in this region are due to nuclear binding and nucleon momentum 
distribution (Fermi motion). At small x the effects beyond the lA should be taken 
into account, such as corrections due to scattering on nuclear pions (meson exchange 
currents) and coherent multiple interactions of intermediate hadron (or quark-gluon) 
states with bound nucleons that cause nuclear (anti) shadowing effect. Summarizing, 
for the nuclear structure function (to be specific we discuss F2) we have different 
contributions: 

F^ = F^^"-^ + 5^F^ + d^o^Ft (1) 

where the first term in the right-hand side stand for the impulse approximation, and 
5:iF2 and 5coh^2 are the corrections due to scattering off nuclear pion (meson) field and 
coherent interaction of the intermediate virtual boson with nuclear target, respectively. 

The nuclear SF in impulse approximation can be written as the sum of the proton 
(t = p) and neutron (t = n) contributions 

F^'^\x,Q^)= £ j[Ak]{\+kjM)^'{e,k)Fi{x\Q\k^) (2) 

where the integration is taken over the four-momentum of the bound nucleon k = 
{M + e, k) and [Ak] = d£d-^^/(2;r)^. The energy and momentum distribution of bound 
nucleons is described by proton (neutron) nuclear spectral function ^P^"^'{e,p), which 
is normalized to the proton (neutron) number in a nucleus. In the integrand FjP^"> are the 
structure functions of bound proton (neutron), which depends on the Bjorken variable 
x' =x/{l + {e+k^) /M) , momentum transfer square Q^ and also on the nucleon invariant 
mass squared k^ = {M + eY — k^. We also note that Eq.(2) is written for the kinematics 
of the Bjorken limit assuming that the momentum transfer is antiparallel to z axis. For 
the derivation and more general expressions valid at finite Q^ beyond the kinematics of 
the Bjorken limit see [7]. 

Although the calculation of Fermi motion and nuclear binding corrections with real- 
istic nuclear spectral functions provide a general trend of observed behavior of the ratios 
M2 = F2/F2 (EMC ratio) at large x [7], the quantitative description of data is missing in 
lA. This indicates that the modification of the nucleon structure functions in nuclear en- 
vironment is an important effect. In Ref.[7] this effect is related to off-shell dependence, 
i.e. the k^ dependence, of structure functions in Eq.(2). Since characteristic energies and 
momenta of bound nucleons are small compared to the nucleon mass, the off-shell effect 
can be treated as a linear correction in the parameter v = {k^ —M^)/M^ to the structure 



function of the on- shell nucleon: 

Fiix.Q'^^) =F2{x,Q')[\+v5f2{x,Q')]. (3) 

The function 5/2 describes the relative off-shell effect. In Ref.[7] this function was 
phenomenologically extracted from analysis of data on the ratios ^2- 

Although the lA with off-shell correction gives the leading contribution to nuclear SF, 
it is by no means complete. It is known that the balance equation for the nuclear light- 
cone momentum is violated in the impulse approximation: because of nuclear binding 
the nucleons do not carry all of the light-cone momentum of the nucleus. This indicates 
the presence of non-nucleon degrees of freedom in nuclei which carry the missing light- 
cone momentum. One possible solution is to consider the scattering off the virtual 
nuclear pions that results in a correction to the structure function, 8nF2 (for a review 
and more references see [10, 11]). In effective ;rN theory with pseudoscalar coupling 
it can be shown that the nuclear pion distribution indeed reconciles with the balance 
of the total nuclear light-cone momentum [12]. In Ref.[7] the nuclear pion correction 
was calculated in terms of the convolution of the light-cone distribution of nuclear pion 
excess in a nucleus and is the pion structure function leading to a positive correction at 
small values oix (about 5% correction at the maximum at x ~ 0. 1 for heavy nuclei). 

At smaller x {x <^0.\) the coherent effects in nuclear DIS are relevant (see, e.g., 
[10]). These effects are associated with the propagation of intermediate hadronic states 
in a nucleus. In particular, the interference of multiple scattering terms and the energy 
dependence of scattering amplitudes can cause both the negative (shadowing) and pos- 
itive (antishadowing) corrections depending on the region of x. The rate of this effect 
depends on the scattering amplitude of the virtual hadronic states off the nucleon in for- 
ward direction. In Ref.[7] the interaction of virtual hadronic states with the nucleon was 
described in terms of phenomenological scattering amplitude aeff • The corresponding 
nuclear amplitude and correction to the structure function 5coh^2 was calculated using 
the Glauber-Gribov multiple scattering theory. 

The outlined approach was applied in the analysis of data on the ratios ^2(^/5) = 
F^ 1^2 *^f structure functions of two nuclei [7]. We perform a fit to experimental data 
by minimizing j} = Y,{M^^ —M^)^/o^{M^^), where o^ represents the uncertainty of 
experimental data and the sum includes all the data points. The analysis includes data on 
M2{-A/B) for a variety of targets from D to Pb with Q^ > \ GeV^ and in a wide region 
of X (overall about 560 data points). In this way we phenomenologically constrain the 
off-shell function 8/2 and the effective amplitude agff from data. To this end we used a 3- 
parameter model to parametrize the x dependence of the off-shell function 8/2 and one 
phenomenological parameter to describe the Q^ dependence of effective cross section 
Ceff = 2Imaeff. This approach leads to a very good agreement with data, reproducing 
the observed x, Q^, and A dependencies of nuclear effects in the structure functions, 
with overall ;f^/d.o.f. = 459/556. In order to test the model we performed a number 
of fits to different sub-sets of nuclei in the region from ^He to ^*^^Pb. The results are 
compatible within the uncertainties with the result of the global fit thus indicating an 
excellent consistency between the model and the data for all nuclei. 



Nuclear corrections in neutrino DIS 

In contrast to the CL case, the neutrino interactions with nuclei involve both the vector 
and the axial current and charge exchange processes. For this reason the results on CL 
structure functions can not be automatically taken over the neutrino case. Nevertheless, 
there is certain similarity in nuclear scattering mechanisms for high-energy CL and 
neutrino interactions. In particular, Eq.(2) can be applied to calculate neutrino-nuclear 
SF in impulse approximation. In the analysis of neutrino and antineutrino SF for heavy 
nuclei it is useful to explicitly separate the isoscalar and isovector parts in Eq.(2) by 
introducing the isoscalar, ^^ + ^" =A^o, and the isovector, ^p - ^" = (Z-N) 0^\ , 
nuclear spectral functions. The Z and A^ are the proton and neutron number and the 
reduced spectral functions <^o and !^\ are normalized to unity. Then we can write Eq.(2) 
as 

7A/4 _ 1 / T7P+n\ P I -nP-n 



where i labels the SF type, Ff " = Ff ± F" and the parameter /3 = {Z — N)/A describes 
the excess of protons over neutrons in a nucleus. In Eq.(4) we use the contracted 
notations of the integration in Eq.(2) with reduced spectral functions <^o and ^\, 
respectively. In the kinematics corresponding to the Bjorken limit the integrands in 
Eqs.(2) and (4) are similar for different types of SF (/ = r,L, 2, 3). However, the finite 
Q^ effects in nuclear convolution depend on the SF type. In addition, the effect of the 
nucleon transverse momentum leads to the mixing of Fj and Fi contributions to order 
p\/Q^ (for more detail see [8]). 

Equation (4) has a convenient form for the calculation of neutrino SF with definite C- 
parity. In particular, we apply it to v ± v combinations. For the C-even combination we 
have i^.(^+^)P = i^v+^)", because of the isospin symmetry. For this reason the isovector 
term in Eq.(4) should vanish. However, one should remark that the isospin relations for 
the structure functions for neutrino CC scattering are violated by the mixing of different 
quark generations and the c-quark mass effect, even in the presence of exact isospin 
symmetry for the parton distributions. This effect results in the nonzero difference 
^.(v+v) (/?-«) ^ g|j^2 Q^ ^jj-j^ Q^ j-j^g Cabbibo mixing angle. Nevertheless, in heavy nuclei 
the contribution of isovector term is accompanied by a small parameter /3 and, therefore, 
it is a good approximation to keep only the isoscalar term in Eq.(4). 

The application of Eq.(4) to the V— V asymmetry of F2 and xFt, requires somewhat 
more attention. We first consider ^2^^^ (similar discussion also applies to Fj and Fi). 
This structure function is C-odd and dominated by the isovector quark distributions. The 
isospin symmetry suggests Fj^^^'P = —FJ^'^". Therefore, the first term in the right 
side of Eq.(4) vanishes and the nuclear structure function is determined by the second 
(isovector) term. However, as mentioned above, the SF isospin relations are violated and 
fjv-v)(p+n) ^ sin^0c. Surprisingly, this effect should not be neglected in the analysis 
of nuclear corrections for the difference v — v even in the first approximation. This is 
because the isovector contribution in Eq.(4) is suppressed by the factor J3 and thus the 
relative contribution ofFJ'^^^'^P^"' is enhanced [8]. The difference xF^^^ is C-even and 
includes the (isovector) contribution from the light quarks and the 5-quark contribution 
to the isoscalar part. For the isoscalar nucleon the difference xFj^^'^'^p^"-' is driven by 



the 5-quark distribution [in the leading approximation in the parton model 2x{s + s)]. 
Because /3 is a small parameter in heavy nuclei, the contribution from the isovector term 
in Eq.(4) is suppressed and, therefore, xF^^^ in nuclei is dominated by the strange quark 
contribution. 

Similar to the CL case, the lA should be corrected for the nuclear pion effect (PI) and 
the effects associated with multiple coherent interactions of intermediate hadronic states 
in nuclear environment (NS). The PI correction can be treated similar to the CL case. 
In the convolution approximation we find that the PI corrections to neutrino and CL 
F2 are similar, and the corresponding correction to xF^ can be neglected. The coherent 
nuclear correction 5coh^ essentially depends on the type of the structure function / and 
on the C-parity [8]. In general, the NS corrections 5coh^2'3^^ can be described in terms of 
four effective amplitudes a^l^\ where / = 0, 1 corresponds to the isoscalar (isovector) 
nucleon configuration and C = ±1 is C-parity. For comparison, the corresponding CL 
correction is determined by a single a^^^'^^ amplitude. 
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FIGURE 1. Different nuclear effects calculated for ^"'Pb at g = 5 GeV^ for C-even and C-odd 
combinations of neutrino SF (per one bound nucleon). Shown are the ratio Fj^^^'^/f!^' ^ (left panel) 
andFj '^/(jSFj ) (right panel). The labels on the curves correspond to effects due to Fermi motion 

and nuclear binding (FMB), off-shell correction (OS), nuclear pion excess (PI) and coherent nuclear 
processes (NS). 

Figure 1 illustrates different nuclear corrections to the C-even and C-odd SF ^2^"^^ of 
^'^^Pb calculated in our model [8]. The result for ^2^^^ resembles the EMC effect for 
Fr^ (for comparison, see Fig. 3 of Ref.[7]). At large x the nuclear correction is driven 
by FMB and OS effects. We recall that the off-shell effect, which modifies the structure 
functions of bound nucleon, is important to change the slope of the EMC ratio bringing 
it close to the CL data. At intermediate x values we observe some cancellation between 
the PI and the NS corrections. This latter is the dominant effect for x < 0.05, providing 
a suppression of the nuclear structure function. 

The corresponding nuclear corrections for the difference Fj^^^ are shown after nor- 
malizing them to /3 = {Z — N)/A. The shape and the magnitude of nuclear corrections 
are somewhat different from the C-even case. Note that the coherent nuclear effects are 
particularly pronounced in the C-odd channel, leading to a larger shadowing effect at 



X < 0.005. The enhancement at 0.02 < .<: < 0.2 (antishadowing) with the maximum at 
X ~ 0.05 is a combined effect due to nuclear smearing (FMB) and a constructive inter- 
ference in coherent nuclear correction which is driven by the real part of the effective 
amplitude a[|^~^ 

Neutrino DIS sum rules for nuclei 

The DIS sum rules reflect the underlying symmetry of interactions and put rigid 
normalization constraints on the structure functions. In the study of nuclear DIS the 
sum rules play a special role providing a link between different nuclear effects. 

In this Section we briefly discuss the Adler and the Gross-Llewellyn-Smith (GLS) 
sum rules for nuclei (for more detail see [8]). The Adler sum rule provides the relation 
between the integrated difference of the isovector combination ^2^ ^ to the isospin of 
the target [13]: 

Sa= dx F^-\x, Q')/{2x) = 2/„ (5) 

where the upper integration limit, the ratio of the nuclear and the nucleon masses M^/M, 
is the maximum possible value of the Bjorken variable x = Q^ /{IMqo) for the nuclear 
target and 4 is the projection of the target isospin vector on the quantization axis (z axis). 
In the quark parton model the Adler sum is the difference between the number of valence 
u and d quarks of the target. The Adler sum rule is valid beyond the parton model and 
survives the strong interaction effects because of the underlying isospin symmetry of 
strong interaction. However, in the derivation of the Adler sum rule the effects of non- 
conservation of the axial current that limits its applicability to the region of relatively 
large Q^. Also Eq.(5) should be corrected for Cabibbo mixing effect (see, e.g., [9]). 

For the proton Sa{p) = 1 and for the neutron 5a(«) = —1. For a generic nucleus of 
Z protons and A^ neutrons the Adler sum rule reads (we consider the nuclear structure 
functions per one nucleon): 

Sa={Z-N)/A = [5. (6) 

We now discuss in turn the contributions to 5a from different nuclear effects. First we 
calculate the Adler sum rule for nuclear SF in impulse approximation. By performing 
the direct integration by x of Eq.(4) and using the isospin relations for the proton 
and neutron structure functions we verify that the FMB correction exactly cancels out 
and we obtain (6). Also, the nuclear pion correction to 5a vanishes because of charge 
conservation. However, the OS and NS contributions to 5a are both individually finite. 
The requirement of cancellation between these corrections provides important constraint 
on the isovector part of nuclear effects. We note that the NS correction to Fj^^ depends 
on effective scattering amplitude a^^^' and a^\f^^. The amplitude a^^^^' also determines 
the NS corrections to F^ and was phenomenologically determined in the analysis of [7]. 
The Adler sum rule helps to fix the Re/Im ratio of a^\}^^ [8]. Note also that the Adler 
and the GLS sum rules were used in [8] in order to evaluate the effective cross section 
aeff = 2Ima(^^+'' governing the NS effect at high Q^. 



The GLS sum rule is the integrated structure function F^^ 



5'gls 



Ma/M 



daF^+^x,Q^)/l, 



(7) 



where we write the GLS integral per one nucleon for a generic nuclear target. In the 
quark parton model the GLS sum gives the number of valence quarks of the target 
"^GLS — 3 [14]. Since the number of valence quarks is directly related to the baryon 
number, in the parton model 5gls is independent of nuclear target. However, in QCD 
the direct relation between the baryon current and 5gls only holds in the leading twist 
approximation (LT) and in the leading order (LO) in strong coupling constant ag. In 
contrast to the Adler sum rule, 5gls depends on ^ and is affected by QCD radiative 
corrections, target mass and the higher-twist effects. In general, 5gls also depends on 
nuclear target. 

We explicitly separate nuclear corrections to the GLS integral as 55gls 
5qls' where Sq^s refers to the GLS integral for the nucleon (proton) and calculate the 
GLS integral (7) for a number of nuclei as a function of Q^ [8]. The results are shown in 
Fig. 2. 
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FIGURE 2. The GLS integral 5gls (upper panel) and its variation with respect to the nucleon 55gls 
(lower panel) calculated in [8] for different nuclear targets as a function of Q^. Data points are CCFR 
extraction of the GLS integral for iron [16]. 

A few comments are in order. It can be verified analytically using Eq.(2) that 55gls = 
in lA [15], i.e. the FMB correction to the GLS integral vanishes, similar to the Adler 
sum rule. The cancellation between phenomenological OS and NS corrections at high 
Q^ was verified in [7] (the normalization of nuclear valence quark distribution was one 



of the constraints of analysis [7]). However, at low Q^ the cancellation between different 
nuclear effects in 5gls is not exact. 

The results of Fig. 2 were calculated with the LT nucleon SF using the NNLO 
coefficient functions and the NNLO parton distributions of [17] and also include the 
target mass correction. As can be seen from Fig. 2 the correction 55gls is negative 
and of the order of 2% at low Q^ and decreases progressively with Q^ and the results for 
different targets are similar. This is mainly due to the explicit (^ dependence of effective 
cross section describing the nuclear shadowing effect. We note that phenomenological 
cross section effectively incorporates contributions from all twists since it is extracted 
from data. Higher twists are known to be important at low Q^ and for this reason we 
should not expect the exact cancellation between the OS and NS corrections in this 
region. We also note that the general Q^ dependence for ^^Fe is in agreement with the 
CCFR measurement [16]. 

The nuclear effects somewhat modify the Q^ evolution of the GLS integral. For this 
reason nuclear corrections should be taken into account for precise extractions of a^ 
from the GLS sum rule. 
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